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Abstract 

Non-equilibrium spin transport through an interacting quantum dot is analyzed. The coherent 
spin oscillations in the dot provide a generating source for spin current. In the interacting regime, 
the Kondo effect is influenced in a significant way by the presence of the precessing magnetic field. 
In particular, when the precession frequency is tuned to resonance between spin up and spin down 
states of the dot, Kondo singularity for each spin splits into a superposition of two resonance peaks. 
The Kondo-type cotunneling contribution is manifested by a large enhancement of the pumped 
spin current in the strong coupling, low temperature regime. 
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Introduction. — The rapid progress of nano-electronics and information technologies has 
prompted intense interest in exploiting the interplay of electron charge and spin properties, 
which results in the emergence of semiconductor spintronicsQ, 0]. Since it is a manageable 
task to manipulate and measure electron and nuclear spins via optical or electric field ^ , 
the solid-state implementation of a quantum information processing architecture is intensely 
studied j^, aiming at the possibility of large-scalability of the quantum computer. One of the 
most important spin-based electronic devices is a mesoscopic quantum dot (QD) system, in 
which the spin co;e.„ce fa elect., o. nucM . .Miv., ,o„ JQ. Sp;.pola..ed 
transport throiigh a QD has been extensively investigated recently. It has been shown 
theoretically 0, 0| and demonstrated experimentally that a QD system will function as 
a phase-coherent spin pump in the presence of sizable Zeeman splitting. A QD-based spin 
battery device has been proposed by use of a combination of a vertical magnetic field and 
an oscillating electric field 10| . Very recently, spin-polarized current has been detected from 

. For the latter, 

due to the large charging energy at low temperatures, a more subtle effect — the creation 
of new states of many-body character at the Fermi level by the Kondo effect — is expected 
to occur 13] . Although the effect of Kondo resonance on the charge current has been well 
studied in a QD|lJ], its influence on the spin current and spin detection might be equally 
significant, since it provides a spin- flip cotunneling mechanism. 

In this paper we explore spin current response to a magnetic-oscillating quantum dot, 
which is schematically shown in Fig. 1. The setup is similar to the one used in recent 
ESR-STM experiments [IS]. The single-electron levels in the dot are split by an external 
magnetic field Bq, eai — £dT: = QI^bBq, where g is the effective electron gyromagnetic factor 
and /iB is the Bohr magneton. The two spin levels are coupled by a rotating magnetic field 
(Sicoscut, Bismut)^ applied perpendicularly to the field Bq. The possibility of the ESR 
setup based on GaAs QD has been discussed in Ref. [9]. The QD is in some sense similar 
to an oscillating magnetic dipole moment of a two-level atom. However, we show that the 
spin transport properties are highly non-trivial for the QD system, due to the many-body 
excitations and coupling with an external electrode. 

Theoretical formalism. — The model Hamiltonian under consideration can be written as 
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H = J2 ^dadtd^ + Ud^^d^d^di - gfiBBi{d+die'''' + d+d^e-'""') (1) 

+ Yl ^k4a(^ka + [^^ataJa + H.C.] . 

ka ka 

Here d'^ and a^^ create an electron of spin a in the QD and in the lead, respectively. The 
first two terms in (1) represent correlated spin levels of the QD, the time-dependent terms 
describe the coupling between two QD spin states | t) and | i) by the rotating magnetic 
field, and the last term is the hybridization of the QD to the lead. A similar model has been 
proposed for electrical detection of single-electron spin resonance jl?!. llS^ . 

The time dependence of the Hamiltonian (1) may be eliminated by introducing a uni- 
tary transformation U = e-^'^E*("fci"'=i-"feT'^ftT)+(4'^i-4'^T)] jg ^et to be unity) and thus 
redefining the Hamiltonian in the rotating reference (RF) as follows 

Hrp = U-'HU + i—^U (2) 
= "Y^^dad^dcj + Ud'^d^d'^di — gfj,BBi{d^di + d'J^di) 

a 

+ Y ^kaat^ttka + J2 [^kat^ada + H.C.] , 
ka ka 

where Eda = £da ± f and eka = Cfc ± f are the dressed dot and lead levels for up and down 
spins, respectively. One can see that the rotating magnetic field, in effect, shifts lead electron 
energy to the opposite directions for up and down spins. This spin-oriented energy splitting 
originates from the hybridization term in (2), which propagates energy from the QD to the 
lead. Thus the spin- up chemical potential is now /i^ = —uj/2 (The reference energy is set 
Ep = 0), while the spin-down chemical potential is fi^ = uj/2. The charge chemical potential 
is still fie = (/^t + = 0. The distinction between /Xo- and /Xe was recently discussed in 

Ref.[19]. 

The physical picture of spin current generation is as follows. Due to the one-photon 
difference between the spin-down and spin-up chemical potentials of the lead, an electron in 
the spin-down channel can tunnel into the Edi level in QD. After a spin- flip process given by 
Bi term in (2), the electron can tunnel out of the QD into the spin- up channel. This process 
repeats, so a steady spin current builds up in the lead. Note that the net charge current is 
zero because there is only one electrode in the system. This zero charge current feature is 



retained when more electrodes are taken into account with no bias voltage applied between 
them. We emphasize that the present situation is not a simple realization of magnetization 
of the electrode by the rotating magnetic field. Here the generation of spin current is due 
to real charge transfer processes between the two spin channels of the lead. Therefore, the 
spin current can be conveniently controlled in experiments. To calculate the spin current 
we make use of the non-equilibrium Green function formalism. 

In the absence of hybridization and Coulomb interaction, the QD is reduced to a simple 
driven two-level system (TLS). Its dynamics is featured by a coherent weight transfer (Rabi 
oscillations) between the two spin states, which is complete when the rotating frequency is 
tuned to the resonant condition uo^i = q^bBq. The spin oscillation period is T = Tr/fi, where 



Vt = -^/A^ -|- A^gjiBBi)"^ is the Rabi frequency and ^ = ujr — uj denotes the detuning from 
the resonance. 

To proceed we introduce a canonical transformation 
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where (f) = tan-^(^f^). In terms of new fermion operators, the dot Hamiltonian in Eq. 
(2) is has a diagonal form '^^ScaC^Ca + Uc^c-^c^ci with — {^di + ^di i ^^)/2 for the up 
and down spins. 

The spin-o", say spin-f, current can be calculated from the time evolution of the occupation 
number N^{t) = X^feC^^Cfci for electrons in the electrode using the non-equilibrium Green 
function: 

^T = -e|(iVTW) = |Rc|j]V^, I ^[uG<,(6)]n|, (4) 

where we have defined the lesser Green function matrix [G^^(t)]o-cr' = i{c^i{0)O'kcrit)) ■ Next, 
we use Dyson's equation to calculate G^^ and express by the Green functions of the dot 
as follows 

= f / ("^^^^^'^ + fMGlie) - G^(6)]}u^)^^ , (5) 

where r(e) = 27r^^|Vfcp5(e — e^) and [G'^{t)]„„i = '^(c^/(0)co-(t)) are the lesser Green 
functions of the dot, A = G^ — G" is the spectral function matrix, and /^(e) — /(e — l^a) 
are the Fermi distribution functions for the spin-up and spin-down channels. 

The non-interacting spin current. — In the non-interacting {U — 0) case or within the 
mean field (MF) treatment of the Coulomb interaction for finite U, the retarded Green 



function matrix is given by 

where Sq = —iTI is the tunnehng contribution to the retarded self-energy, and g,l{e) is the 
unperturbed Green function matrix of the dot, which in the mean-field approximation has 
a simple Hartree-Fock type [g^{e)]aa' = ^aa'ijz^^zu + with = (c+Co-). must 

be calculated self-consistently via the relation Uar = Im J ^[^^(e)]^-^- The lesser Green 
functions in (5) are obtained via the Keldysh equation — G^Sq G" with the lesser 
self-energy matrix 

Y,< =iT\ ^^^^ ^^^^ ^^"^ ^ ^™ ^ ^^^^ ~ ^^^^ 1 (7) 

y sin0cos0[/|(e) -/t(e)] /|(e) sin^ + /^(e) cos^ j 

Substituting the expressions for G^ into (5), we obtain the expression for the tunneling spin 
current 

"^^ ^ f / ^{"GS(e)S<GS(6)u^}n{/T(6) - A(6)}, (8) 

where 

< / sin^ 6 sin 6 cos 6 \ 

So = r ! • 

\ sin (j) cos (j) cos (j) j 
Since the electrons are flowing from spin-down channel to spin-up channel, we define the 
total spin current as Jg — Ji — J^; — — 2J|. Although the coherent spin superposition is 
strongest at two-level resonance ui = cor, the spin current is also influenced by the spin 
levels £do-relative to /ig-. Figure 2(a) shows the non- interacting spin current versus rotating 
frequency at different temperatures. The undressed spin levels are set to lie deeply below 
Ep. One can see that at low temperature (solid line), the spin current peak is not at, but far 
from u! — ujR. This is because at a; = ujr the dressed spin- up level is still lower than /x-f 
and thus the electron in the dot cannot tunnel out via small thermal excitation. Due to the 
state exchange in QD, as shown in the inset of Fig. 2(a), when to increases further, crossing 
the avoided crossing between the two spin levels, the upper eigenstate is dominated by spin 
component | |) while the lower state by | [). Thus when the upper level is higher than 
spin tunneling processes occurs, and a peak develops in the current spectrum. When the 
temperature is increased to ksT ~ — the spin-current peak shifts io uj — ujr with a 
larger peak amplitude. Further increasing temperature such that ksT > cu will smear out 
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the mismatch between yU| and and thus the spin current amphtude begins to decrease, 
as shown in Fig. 2(a). The dot levels can be conveniently controlled by the gate voltage Vg. 
When Vg is modulated such that are at the interval between fj,^ and then the spin 
current peak will occur exactly at = uji [see Fig. 2(b)]. 

The model can also illustrate the spin-current generation in the Coulomb blockade regime 
when Coulomb interaction U is treated in the mean-field formalism. As an example, we show 
in Fig. 3 the spin current versus gate voltage for several values of uj. Two Coulomb peaks 
with interval U can be resolved, and as in Fig. 2, the current amplitude increases when the 
magnetic field frequency is tuned towards to two-level resonance u = ujr. 

The electronic correlation effect. — In the low temperature, strong coupling, and large 
charging energy regime, correlation effects enter and are expected to significantly infiuence 
the spin transport. To illustrate correlation effects, we again use the equation of motion 
method to solve the retarded Green function in (2), which generates higher-order Green 
functions that have to be truncated to close the equation^^l- After a straightforward cal- 
culation, in the infinite-[/ limit we obtain 

G'-(e) = [el-% - SS(e) - ^\{e)]-\l - n,), (10) 

where isdjaa' = ^aa'^da a^d the interaction self-energies are given by [S^(e)]o-o-' = 

5<r^'Y.k[¥r^T^ + rl^'^'f^+r^ ]- (nc).<x' = 5..'(4c^) needs to be calculated self- 
consistently via the relation (c^TCo-) = Im J ^G^^(e). G'^(e) is difficult to obtain since 
the lesser self-energy cannot be given exactly in the interacting regime. However, the sit- 
uation simplifies in the steady-state transport [2I | , in which the occupation on the dot is 
time-independent, i-^{c'^{t)ca{t)) = 0. Writing down the Heisenberg equations for c^{t) and 
Co-(t), and integrating out the lead electron operators, we find the self-consistent integral 
equations for the retarded Green functions as follows: 

(c+ct) = Im^ ^G< (6) = J^ Im(G^^ - Gy^)[cos2 0/t + sin^ (11a) 

(c+q) = Im J ^G<(e) = J ^lm{Gl, - G[,)[sin' <Pf^{e) + cos' <Pf,{e)], (lib) 

which now, together with (10), close the equations of motion for G*^. Similarly, the quantities 
(c^Cct) = —i J ^G'^g:(e) can also be expressed by the retarded Green functions. Thus we 
can calculate the spin current in (5) directly after G^'(e) is numerically obtained. 
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The magnetic-field-induced excitation properties are shown in Fig. 4 by plotting local 
spin-resolved spectral densities Pa{€) — [uG^(e)u'^]o-(T for different values of a;. In the absence 
of a rotating magnetic field {uJ,Bi = 0), as observed from Fig. 4(a), the spectral density 
for each spin component is characterized by a broad single-particle peak around Eda, and 
a sharp Kondo peak at Zeeman energies e = —oon for up spin and e = ujr for down spin. 
The spectral weight of p| is enhanced by the magnetic field, while pi is greatly suppressed. 
Therefore, a net spin moment develops in the QD. Physically, the Kondo peak at e = —ujr 
is due to the cotunneling processes (see inset) in which one spin-down electron in the lead at 
the Fermi level tunnels into the dot and occupies energy level £|, followed by another electron 
on e-i tunneling out to the spin-up channel with the energy ujji below from Fermi level Ep. 
Such cotunneling processes transfer electrons from spin-up channel to spin-down channel 
of the lead. Thus a spin current may be induced in an efficient way by these Kondo-type 
cotunneling processes. Similarly, the Kondo peak at e = ujr in Pj(e) is due to a superposition 
of many cotunneling processes in which one spin down electron with the energy ujr above 
the Fermi level tunnels into the dot and simultaneously another spin-up electron in the dot 
tunnels out to the Fermi level of the lead. Both kinds of cotunneling processes contribute 
to the spin-current generation. When the rotating magnetic field is switched on, the spin 
Kondo resonances and spectral structures begin to change. At a; < a;ij, the broad single-spin 
and sharp Kondo peaks shift toward higher energies for up spin, while they shift to lower 
energies for down spin [Fig. 4(b)]. At two-level resonance u — ur, as shown in Fig. 4(c): (i) 
each spin Kondo singularity splits into two prominent peaks around /i^ and //j; (ii) P|(e) and 
Pl(e) completely overlap with identical spectral weight. This overlap is caused by the fact 
that dXoj = ujr, each spin eigenstate in the dot is a strong superposition of | ]) and | |) states 
with equal probabihty amplitude. When uj increases further, the two-level superposition is 
suppressed, and the state exchange occurs. As a response, the spectral densities of two spins 
also exchange their structures; as can be seen from Fig. 4(d), pj.(e) is enhanced while p^{e) 
is greatly suppressed. 

To see the influence of the many-body correlation effect on the spin transport, we show 
in Figure 5 the spin current as a function of rotating frequency. Compared to the non- 
interacting case, it can be seen that the spin current is greatly enhanced by the cotunneling 
processes in the strong coupling, low temperature regime. A detailed study of cotunneling 
effects on the spin current is in progress. 
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In summary, we have analyzed the spin current properties of an interacting QD exposed 
to a rotating magnetic field. The spin-fiip process in the QD and the effective difference of 
two spin chemical potentials suffices the generation of spin current with no charge current. 
The interplay of the many-body cotunneling process and coherent two-level resonance in the 
QD produces a remarkable enhancement of spin current. Since a controllable generation of 
spin current provides an efficient source for spin injection, spin detection, and virous variety 
of spin-based devices, we expect the present results may have practical applications in the 
field of spintronics. 

This work was supported by grants CNSF 90103024, MOST-G001CB3095, and U.S.- 
DOE-DE-FG03-98ER45687. 
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Figure captions 

Fig. 1. Model setup for spin transport through a QD. The degenerate spin levels in the 
dot are split and coupled by the external magnetic fields Bq and Bi, respectively. Due to 
the imbalance between spin-up and spin-down chemical potentials (see text), the electron 
tunnels from the spin-down channel, through the QD, to the spin-up channel. 

Fig. 2. (a) Non- interacting spin current as a function of rotating frequency at different 
temperatures with the undressed spin levels set = —2.5, e^i = —1-5; (b) Contour plot 
of spin current as a function of gate voltage and rotating frequency. The spin levels in the 
dot are modulated by the gate voltage as EdaiVg) = Sda + eVg with Ed] = 1, Sdi = 2. Other 
parameters are F = 0.1 and q^bBi = 0.2. 

Fig. 3. Spin current versus gate voltage for different rotating frequencies with U = 2. 
Other parameters are the same as in Fig. 2(b). 

Fig. 4. The spin up (solid line) and spin down (dashed line) spectral densities at different 
rotating frequency. Parameters are Ed^ = — 2.5r, Sdi = — l.SF, T = O.OOIF, and qubBi = 
0.2T. 

Fig. 5. Spin current as a function of rotating frequency in the Kondo-type cotunnel- 
ing regime (solid line). The non- interacting spin current is also shown (dotted line) for 
comparison. Other parameters are the same as in Fig. 4. 
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Fig. 1 Zhang et al. 




Fig. 2 Zhang et al. 




Fig.4 Zhang et al. 
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Fig. 5 Zhang et al. 



